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Abstract 

The propagation of the scalar (a) and vector (a;) mesons in nuclear matter 
is studied in detail using the Walecka model over a wide range of densities and 
including the effects of a finite a width through the inclusion of a two-pion loop. 
We calculate the dispersion relation and spectral functions of the a and (transverse 
and longitudinal) u mesons, including the effect of a-to mixing in matter. It is shown 
that the mixing effect is quite important in the propagation of the (longitudinal) u 
and a mesons above normal nuclear matter density. We find that there is a two-peak 
structure in the spectral function of the a channel, caused by a-uj mixing. 

PACS numbers: 21.65.+f, 24.10.Jv, 21.30.Fe, 25.75.-q, 14.40.-n 

Keywords: quantum hadrodynamics, dense nuclear matter, vector-scalar mixing, disper- 
sion relation, spectral function 
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Recently much attention has been focussed on the variation of hadron properties in hot 
and/or dense nuclear matter]!]. In particular, the medium modification of the light vector 
(p, D and 0) meson masses has been investigated by many authors@, [|, |J ^ ||, 0, |], |§. 
Experimentally, dilepton production in hot and/or dense nuclear matter produced by 
relativistic heavy ion collisions can provide a unique tool to measure such modifications. 
It is well known that data obtained at the CERN/SPS by the CERES Q and HELIOS 
collaborations show a significant amount of strength below the free p meson peak. Some 
authors [O have concluded that this is caused by a downward shift of the rho meson 
mass in dense nuclear matter. To test this idea, measurements of the dilepton spectrum 



from vector mesons produced in nuclei are planned at TJNAF[13|] and GSI||14|| (see also 
Rat 

On the other hand, the a meson has been treated as a correlated two-pion state in 
the scalar channel for a long time. However, recently some people have reanalysed the 
s-wave (1=0) 7T7r-scattering phase shift, and argued for the existence of the a as a genuine 
resonance [HJ . From the theoretical side, it is also anticipated that the a meson may play 
an important role as the chiral partner of the 7r meson |T^| . Some experimental possibilities 
have been proposed for investigating the behaviour of the a in hot and/or dense nuclear 
matter [0, 0. 

It is well known that in nuclear matter the a can couple to the longitudinal mode of 
the uj meson [[0J. This work is an extension and elaboration of studies in Ref.|| where the 
invariant mass of the u meson moving in nuclear matter was studied. Here, we generalize 
this work by including free space widths for the a and uj and by studying a-u mixing 
in some detail as a function of density. Recently Wolf et a/.[[2(| studied the effect of <j-uj 
mixing on the e + e~-pair production in relativistic heavy ion collisions. However they 
calculated only the lowest order mixing diagram, and did not include the nucleon-loop 
diagrams for the a and uj meson propagators. The interesting results obtained in these 
earlier papers and the experimental relevance of the problem suggests that it is time for 
a more complete self-consistent calculation. 
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A complete investigation of the propagation of a and uj mesons in a dense nuclear 
medium requires that one includes the effect of g-uj mixing. In this letter we will study in 
detail the medium modification of o and uj mesons in dense (symmetric) nuclear matter 
using the simplest version of quantum hadrodynamics (QHD or the Walecka model) |2lt| 
and including a-uj mixing, oixix coupling, and the uj width. 

The lagrangian density of the Walecka model (QHD-I) is well known and a full descrip- 
tion of it can be found in Ref.f2]J. We first treat the nuclear ground state in relativistic 
Hartree approximation (RHA), which is also explained in Ref . [ TL\ . Here we shall ig- 



nore the coupling to the channels involving iso-vector mesons, neglecting, for example, 
uj + N-+7r + N&nduj + N^p + N[^. 

To compute the meson propagators, we sum over the ring diagrams, which consist 
of repeated insertions of the lowest order one-loop proper polarization part. This is the 
relativistic, random phase approximation (RPA) [19| . Since we want to include <j-uj mixing, 
it is convenient to use a meson propagator in the form of a 5 x 5 matrix with indices 
a, b running from to 4, where 4 is for the a channel and ~ 3 are for the uj. 

Dyson's equation for the full propagator, T>, is given in matrix form as: 

V = V° + V°UV, (1) 
where T>° is the lowest order meson propagator given by a block-diagonal form as 



D% 



(2) 



\ A 

In Eq.fl^) the noninteracting propagators for the a and uj are given respectively by 

A„(<7) = i \— ~, (3) 

% - m * + k 

DIM = 2 ^ 2 T • , (4) 

where = —g^ u + q 2 ^) , q 2 , = q 2 — \q\ 2 , and m a and ra w are the free a and uj meson 

masses. 
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The polarization insertion in Eq.(|T]) is also then given by a 5 x 5 matrix, 

(5) 



n 



\ II™ (q) U s (q) 

where the lowest order scalar, vector and scalar-vector-mixed polarization insertions are 
given respectively as 

H.(«) = -^/^Tr^ (6) 

IW?) = -igl J -0- 4 Tr[G(k) ltl G(k + q)^}, (7) 

n?(«) = i9&f^T±[G{k) 7ll G{k + q)]. (8) 

Note that we have added the contribution from the pion-loop to II S (the second term in 
the r.h.s. of Eq.(^D) in order to treat the a more realistically. The pion propagator, A w , 
is given by equation (§) with the pion mass m n instead of m a . Here g v , g s and g a7T are 
respectively the nucleon-u;, nucleon-a and a-ir coupling constants. We denote G(k) as the 
self-consistent RHA nucleon propagator, which is given by the sum of the Feynman (F) 
part and the density- dependent (D) part as 



G(k) = G F (k)+G D (k), 



(^k; + M*) 



kf - M* 2 + ie E k 



77T — » 

+ ^-AK-El)0{k F -\k\) 



(9) 



where k*^ = (k° - g v V°, k) (V° is the mean value of the u field), E* k = \J k 2 + M* 2 (M* 
is the effective nucleon mass in matter) and kp is the Fermi momentum. Therefore, each 
polarization insertion (except the pion-loop contribution) can be divided into two pieces: 
the Feynman (F) (or vacuum) part, which does not involve the 9(kp — \k\) term, and the 
density-dependent (D) part. The explicit expressions for the various components of D can 
be found in Ref . |23| . 



We can work out the F parts using the method of dimensional regularization to remove 
the divergence in the loop calculations. We show the results explicitly. For the a, we must 
renormalize two terms: the nucleon-loop and pion-loop diagrams. For the nucleon-loop 



contribution to the a we introduce the usual counter terms to the lagrangianpT 



1=2 



'(J 



(10) 



which includes quadratic, cubic, quartic and wavefunction renormalization. To get the 
"physical" properties of the a meson in free space, we impose the following condition on 



the F part of the nucleon-loop diagram |[24|: 



nf -^(gj, M* = M) = Anf-^(gJ, M* = M) = at q\ 



where M is the free nucleon mass. We adjust the coefficients a 2 and ( to satisfy this 



condition. For «3 and «4 we use the usual values given in Ref . [gjj . For the pion-loop, 
we also introduce an appropriate counter term to the lagrangian, and require a similar 
condition to Eq. (|TTD for the polarization insertion of the pion-loop diagram in free space: 



^eUr loop (ql) = at q 2 = m 2 a . 



'12) 



Finally, we find 

njr-tooP(g) 



2tt 2 



M 



*2 



6 



2 In 



M* 



+ f(x q ) - /(*.)) 



Q 2 ( M* 2 M 2 N 

+ T ^(/K)-2)-^(/(z s )-2) 



9„ 

r*2 /\#2 



(M - M 2 )(f(z s ) - 2) + 2M(M* — M) + 3(M* - M)' 



where i„ 



4M' 2 

„2 , 



^and 



3f7 2 



where = 1 — -^r 1 and 2^ = 1 — ^f-, and 



4m; 



(13) 



(14) 



f(y) 



V^ m T 



?7T 



-7/ tan 



-l l 



for 1 < y < +oo 



y/y, for < y < 1 



for ?/ < 



(15) 



For the u), we add the counter term for the wavefunction renormalization to the la- 
grangian and require the following condition for the F part of the polarization insertion 
of the nucleon-loop: 



K~ l °° P ^ = ^ N - l °°v(q% M* = M) = at q\ = m 2 . (16) 



Then, we get 



W 67r 2^ 



g 2 ,r , M* ( M* 2 M 2 \ , . 

" ~ 2 2 In — -4 — 17 

M \q 2 m l 



11 U) / 

M* 2 \ / M 2 \ 



where z v = 1 — ^M-. Note that vacuum fluctuations do not contribute to the mixed part 
of Eq.(|). 



Using these polarization insertions, we define the dielectric function e as [[03, pi], [23 



det(l -£> II), 



4 x es-L, 



where €t[sl\ is the dielectric function for the transverse (T) [scalar and longitudinal (SL)] 
mode. (Note that the full propagator Eq.([l|) can be rewritten as V = V°/(l - V°Tl) and 
hence detP = detX>°/e.) Then, we find 

e T = l-d U T , (19) 

e SL = (l-d U L )(l-A U s )-iA d (U^) 2 . (20) 

q l 

Here q = \q\ and 1 = q 2 — m 2 + im^T^ , where we add the width of the lo in free space 
(r° = 9.8 MeV). The transverse (IFr) and longitudinal (11/,) components of the polariza- 
tion insertion Eq.(^) are respectively defined by rin (or II22) and II33 — IIoo (we choose 
the direction of q as the z-axis). Note that II™ is the 0-th component of the mixed po- 
larization insertion Eq.(|S]), which vanishes when q = (because of current conservation). 
The eigencondition for determining the collective excitation spectrum is just equivalent to 
searching for the zeros of the dielectric functions. Since we are interested in the medium 



modification of the meson propagation, we restrict ourselves here to the meson branch in 
the time-like region. 

To study the dispersion relation of the meson branch, we first have to solve the nuclear 
ground state within RHA. Because we choose cti and ( to satisfy the renormalization 
condition for the a at q 2 = m 2 a (see Eq. ([TT|) ), the total energy density is written as 



(21) 



4M 2 2" 

where Sq is the usual one (in RHA) given in Ref.[^l|]. (Note that in Ref.||21|| the renormal- 



ization condition on nucleon loops is imposed at <2^=0. This difference yields the second 
term of the r.h.s. of Eq. (^l|) |24| . As a result, our model gives the same physical quantities 



as those of Ref.plj.) Requiring the saturation condition, £/p B — M = —15.7 MeV at 
Po = 0.15 fm -3 (where ps and po are respectively the density of nuclear matter and the 
saturation density), we determine the coupling constants g 2 s and g 2 : g 2 = 101.50 and 
gl = 72.296. In the calculation we fix the values of the free hadron masses to be M = 939 
MeV, m a = 550 MeV, m w = 783 MeV and m w = 138 MeV. This yields the effective 
nucleon mass M*/M = 0.73 at po and the incompressibility K = 452 MeV. We do not 
consider the possibility of medium modification of the pion loop or any of the coupling 
constants in the present work. 

To determine the coupling constant g a7T , we adjust it to fix the width of the a in free 
space. Since the full propagator of a in free space is given in terms of the pion-loop 
polarization insertion, the width T° is given by 

T° = -^ a tql = ml. (22) 

m a 

If we choose T° = 300 MeVgCj, we find g aw = 18.33. 

Now we are in a position to show our results. First, in Fig. [l], we show the dispersion 
relation for the meson branches, which are calculated by searching for the zeros of the 
real part of the dielectric function given in Eq.(^). As illustrated in the figure, if the 
a-uj mixing is ignored, the longitudinal (L) mode and the a (S) mode cross each other at 
one point. However, once the mixing is involved, the two modes never cross each other. 
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Figure 1: Dispersion relation for the scalar and longitudinal modes at Pb/po = 2. The 
curves specified by u are for the longitudinal component of the u meson. The dotted 
curves are the results without a-u mixing. The solid curves include the effect of mixing. 
(Note that the transverse mode is very close to the upper solid curve.) 



We can understand this phenomenon as a level-level repulsion due to the mixing, which 
is quite familiar in conventional nuclear physics, for example, in the Nilsson diagram (see 



Ref.[25Q). Adding the off-diagonal (or mixing) matrix element in the Hamiltonian leads 
to new eigenstates which never cross each other. The new states approach the original 
crossing modes as the off-diagonal part becomes weaker. This occurs at mid and high 
density (above ~ po) because its origin is the mixing effect, which vanishes at ps = 0. 



In Figs. and [| we show the "invariant mass" (m* nv = \Jql — q 2 , with q chosen as 
in Fig. 1 so that the real part of the dielectric function vanishes for that value of q) as a 
function of density. As shown in the figures, we can clearly see the role of mixing, which 
is quite important in determining the meson mass at mid and high density. If the mixing 
is ignored, the masses of the L and S modes cross each other, like the dispersion relation 
shown in Fig. [l|, and the L mass is below the other two at high density. However, in the 
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Figure 2: "Invariant mass" without o-lo mixing. The solid curves are for three-momentum 
transfer q=l MeV. The L and T modes of the u are almost degenerate. The dotted curves 
are for g=500 MeV, in which case the L and T modes are well separated. 



case where the mixing is included the L and S masses never cross each other. At high 
density the L mass is pushed upwards, while the S mass is pulled downwards because of 
mixing (through the level- level repulsion). From the figures we can see that the mixing 
effect becomes vital above ~ p . 

Next we calculate the spectral functions of the T, L and S modes. It is very interesting 
to study these functions because in a thermal model the dilepton yields in heavy ion 
collisions would be proportional to the spectral functions. The spectral function is usually 
defined in terms of the imaginary part of the full propagator. For the T mode it is: 



ST(m* nv ,q,p B ) = ~9fm 

71 



d 



(23) 



1 - d U T _ 

While, in practice, it may be difficult or even impossible to separate the L and S modes 
in an experiment, it is of considerable theoretical interest to study the following spectral 
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Figure 3: "Invariant mass" including a-u mixing. The solid curves are for q—1 MeV. The 
L and T modes of the u are again very close to each other. The dotted curves are for 
g=500 MeV. 



functions: 

S L(m* nv ,q,p B ) = --3m 

7T 

Ss{m* nv ,q,p B ) = --Qm 

7T 

corresponding to the complete, diagonal, longitudinal u and a propagators, respectively. 

In Figs. (|, |^ and ||, we show the shape of the spectral function Si (i=T, L or S), as 
a function of the invariant mass and three-momentum transfer, at ps/pn=2. For the T 
mode the shape is very simple, while for the L and S modes they are complicated because 
of the effect of mixing. In particular, the S mode is quite remarkable: at q=0 MeV there is 
only one peak, while at fixed, finite q there exist two peaks in the spectral function. Note 
that at q=0 MeV the mixed polarization insertion vanishes, so this two-peak structure is 
clearly associated with a-u mixing. As the density grows the two peaks becomes separated 
more widely. In the L mode there is no such structure at fixed q. However, the peak at 
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Figure 4: Spectral function for the transverse u at ps/po — 2. 



low q and that at high q are separated very well and the values of their invariant masses 
are quite different. As the density goes up this tendency becomes more clear. 

Finally, we should comment on the width of the a in free space. In our calculation 
we have used r°=300 MeV|16[ (see Eq.(^)), which is somewhat smaller than usual. We 
have therefore also calculated the spectral functions using a width of 600 MeV but, as the 
results were not qualitatively different we do not report them here. 

In summary, we have studied the propagation of a and u mesons and the effect of 
o-uj mixing in dense (symmetric) nuclear matter, within the Walecka model. We have 
illustrated that the effect of mixing is quite important in the propagation of the longi- 
tudinal uj and a mesons at a density above ~ po- in the scalar (or a) channel we have 
found a two-peak structure in the spectral function at finite three-momentum transfer, 
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Figure 5: Spectral function for the longitudinal uj at Pb/Po — 2. 



which could be measured in future experiments [IB]. It would clearly be very interesting to 
compare these results with the medium modification of the meson propagation using the 
quark- meson coupling (QMC) model [0], in which the effect of hadron internal structure 
is involved. We will report such a study in the near future. 



This work was supported by the Australian Research Council, and the Japan Society 
for the Promotion of Science. K.S. thanks T. Kunihiro, H. Shimizu and T. Hatsuda for 
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Figure 6: Spectral function for the a at Pb/po 
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